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Abstract 



We analyze some exact and approximate solutions to nonlinear equations for heat 
transfer models. We prove that recent results derived from a method based on Lie 
algebras are either trivial or wrong. We test a simple analytical expression based on 
the hypervirial theorem and also discuss earlier perturbation results. 



1 Introduction 

In a recent paper Moitsheki et alp] argued that a method based on Lie al- 
gebras is suitable for obtaining the solution to nonlinear ordinary differential 
equations that appear in simple models for heat transfer. They compared the 
analytical solutions with other results coming from perturbation approaches 
like homotopy perturbation method (HPM) and homotopy analysis method 



tween the users of those fashionable perturbation approaches that arose some 
time ago [6117] . 
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. It is worth noticing that there is an unending controversy be- 
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The purpose of this paper is to determine the usefulness of the results for 
the heat transfer systems provided by the Lie algebraic method and those 
perturbation approaches. In Sec. [2] we analyze the exact solutions arising from 
Lie algebras, in Sec. [3] we outline the application of the well known Taylor- 
series approach, in Sec.HJwe derive a simple accurate analytical expressions for 
one of the models and in Sec. Owe summarize our results and draw conclusions. 

2 Exact solutions 

The first example is the nonlinear ordinary differential equation [1] 



where the prime denotes differentiation with respect to the variable x. This 
equation is trivial if one rewrites it in the following way [(1 + eu)u'\' = 0[1] 
and the solution is 



Moitsheki et alp] derived exactly this result by means of a rather lengthy 
algebraic procedure. It is clear that in this case the Lie algebraic method gives 
us the same answer that we can obtain in a simpler way. 

For the second example 



[1 + eu(x)]u"(x) + eu'{xf = 
u(0) = 1, u(l) = 



(1) 




(2) 



u"(x) -ew(x) 4 = 
u'(0) = 0, u(l) = l 



(3) 



the authors derived the simple analytical expression [TJ 
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They argued correctly that it satisfies u(l) = 1 but they were wrong when 
they stated that "However, u'(0) = only if e = 10/9". Notice that the 
function u(x) = x _2//3 that comes from such value of e does not have the 
correct behaviour at x — 0. Therefore, in this case the Lie algebraic approach 
led to a wrong result. 

Other authors have applied HPM and HAM to the equation [2|3] 

[1 + eu(x)]u\x) + u(x) = 

«(0) = 1 (5) 

with the trivial solution 

hiu(x) + e[u(x) — 1] + x = (6) 

In the following two sections we discuss some of these problems from different 
points of view. 

3 Taylor series 

If the variable of the nonlinear equation is restricted to a finite interval, one 
can try a straightforward power-series solution u(x) = UQ+U1X+U2X 2 + . . . and 
obtain the unknown model parameter from the boundary conditions. In the 
case of the example d2J) the radius of convergence of this series is (e+ l) 2 /[e(e + 
2)] and therefore the approach will be useful for small and moderate values of 
e. As e increases the rate of convergence of the Taylor-series method decreases 
because the radius of convergence approaches unity from above. However, 
this example is trivial and of no interest whatsoever for the application of a 
numerical or analytical method. This reasoning also applies to example (131) 
although in this case we do not have an explicit solution u(x) but x(u). 

The example ([3]) is more interesting because there appears to be no exact 
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solution, and for this reason we discuss it here. The unknown parameter is 
it(0) = u and the partial sums for the Taylor series about x = 

N 

^W = E%KK (7) 

3=0 

enable one to obtain increasingly accurate estimates N increases. Such 

estimates are roots of u^(l) = 1. Although the rate of convergence decreases 
as e increases it is sufficiently great for most practical purposes. Notice that 
the HAM perturbation corrections for this model are polynomial functions of 
x[2] whereas the HPM has given polynomial functions of either x[3\ or e _x [l]. 
However, there is no doubt that the straightforward power-series approach is 
simpler and does not require fiddling with adjustable parameters [2|5] . 

4 The hypervirial theorem 

The analysis of the nontrivial equations for heat transfer models may be easier 
if we have simple approximate analytical solutions instead of accurate numer- 
ical results or cumbersome perturbation expressions. In the case of the models 
(CQ) and (j3J) there is no doubt that the exact analytical expressions should be 
preferred. For that reason, in what follows we concentrate on the seemingly 
nontrivial model ([3]). 

We have recently shown that the well known virial theorem may provide sim- 
ple analytical solutions for some nonlinear problems |8|9] . In particular, we 
mention the analysis of a bifurcation problem that appears in simple models 
for combustion [8]. The only nontrivial problem outlined above is a particular 
case of nonlinear ordinary differential equations of the form 



u"{x)=f(u(x)) 



0<x < 1 
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The hypervirial theorem is a generalization of the virial one. If w(u) is an 
arbitrary differentiable weight function, the hypervirial theorem provides the 
following suitable expression for our problem flS]): 



[w(u)u']'dx = iu(u(l))u'(l) - w(u(0))u'(0) 



dw ^ 
du 



u'y + w(u)f(u) 



dx 



(9) 



In the particular case of the example ([3]) we have 
i 



wiDu'd) 



^V) 2 + ew{u)u A 
du 



dx 



(10) 



When w(u) = u we obtain the virial theorem. Here we also consider the even 
simpler choice w(u) = 1 that we will call hypervirial although it is just a 
particular case. 

Since u"(x) > we try the ansatz 
cosh(6a;) 



cosh(fe) 



that satisfies the boundary conditions in equation ([3]). It follows from equation 
f TTUl) that the adjustable parameter b is a root of 



3e 10fe (56 2 - 2e) + 5e 8fe (126 3 + 9b 2 - lOe) + 30e 6fe (66 3 + b 2 - lOe) 
+30e 4b (66 3 - b 2 + lOe) + 5e 2b (126 3 - 9b 2 + lOe) 
-3(56 2 - 2e) = 

when w(u) = u and 



(12) 



3e 10b (56 2 - e) + 5e 86 (96 2 - 5e) + 30e 66 (6 2 - 5e) + 30e 46 (5e - b 2 ) 
+5e 2fe (5e - 9b 2 ) - 3(56 2 - e) = 



(13) 



when w(u) = 1. 
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Fig. [T] shows u app (0) for some values of e and also the accurate result obtained 
from the Taylor series discussed in Sec. [3j We appreciate that the accuracy of 
the analytical expression fill I) decreases as e increases. However, if one takes 
into account the simplicity of equation (ITT1) the agreement is remarkable. Be- 
sides, the hypervirial theorem with w = 1 proves to be more accurate than 
the virial theorem. It is curious that there is no such test for the HPM or 
HAM 03]. 

As a particular example we consider e = 0.7 (the preferred parameter value for 
both HAM and HPM calculations [2 t 3j). From the partial sums of the Taylor- 
series with N < 30 we obtain m = 0.8186424785. The analytical function fill I) 
yields b « 0.70, u app {0) « 0.80 for w = u and b = 0.657, u app {0) « 0.817 
for w — 1 that is a reasonable estimate of the unknown parameter. Again 
we see that the hypervirial approach is better than the virial one. Fig. [2] 
shows accurate values of u(x) given by the Taylor series with iV = 30, our 
approximate analytical virial expression u app (x) and equation (J3J) for < x < 
1. It seems that the accuracy of u app (x) is somewhat between the HAM results 
of 5th and 10th order [2]. On the other hand, the equation (J3]) derived by the 
Lie algebraic method [1] exhibits a wrong behaviour. 

Finally, in Fig. [3] we compare the numerical, virial (w = u) and hypervirial 
(w = 1) approaches to the function u(x) in a wider scale. We conclude that 
the virial theorem is not always the best choice for obtaining approximate 
solutions to nonlinear problems. 

5 Conclusions 

The purpose of this paper has been the discussion of some recent results for 
the nonlinear equations arising in heat transfer phenomena. The oversimplified 
models considered here may probably be of no utility in actual physical or 
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engineering applications. Notice that the authors did not show any sound 
application of those models and the only reference is a pedagogical article 
cited by Rajabi et al[l]. However, it has not been our purpose to discuss this 
issue but the validity of the methods for obtaining exact and approximate 
solutions to simple nonlinear equations. 

It seems that the particular application of the Lie algebraic method by Moit- 
sheki et al[T] has only produced the exact result of a trivial equation and 
a wrong result for a nontrivial one. Therefore, we believe that the authors 
failed to prove the utility of the technique and it is not surprising that they 
concluded that their results did not agree with the HAM ones [2] (see Fig. [2]). 

We have also shown that under certain conditions the well known straightfor- 
ward Taylor-series method is suitable for the accurate treatment of such non- 
trivial equations. It is simpler than both HAM and HPM [2|3] and as accurate 
as the numerical integration routine built in a computer algebra system[3]. 

Finally, we have shown that the well known hypervirial theorem may provide 
simple analytical expressions that are sufficiently accurate for a successful 
analysis of some of those simple models for heat transfer systems. It is sur- 
prising that our results suggest that the virial theorem [8|9] may not be the 
best choice. 
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Fig. 2. Accurate values of u{x) for e = 0.7 (circles) and approximate results from 
equations (fTTj) (solid line) and (JH) (dashed line). 
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Fig. 3. Accurate values of u(x) for e = 0.7 (circles), virial (solid line) and hypervirial 
(dashed line) approximate results from equation (fTT]) . 
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